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In a half-filled Hubbard model on a square lattice, the next-nearest-neighbor hopping causes 
spin frustration, and the collinear antiferromagnetic (CAF) state appears as the ground state with 
suitable parameters. We find that there is a metal-insulator transition in the CAF state at a critical 
on-site repulsion. When the repulsion is small, the CAF state is metallic, and a van Hove singularity 
can be close to the Fermi surface, resulting in either a kink or a discontinuity in the magnetic moment. 
When the on-site repulsion is large, the CAF state is a Mott insulator. A first-order transition from 
the CAF phase to the antiferromagnetic phase and a second-order phase transition from the CAF 
phase to the paramagnetic phase are obtained in the phase diagram at zero temperature. 



I. INTRODUCTION 



In a hall-filled Hubbard modeli with only nearest- 
neighbor (NN) hopping, the spin exchange between NN 
electrons is antiferromagnetic, leading to the antiferro- 
magnetic insulating ground state, the Neel state, on a 
square lattice. When the next-nearest-neighbor (NNN) 
hopping is present, there is also antiferromagnetic spin 
exchange between NNN electrons, which is destructive to 
the antiferromagnetic order. This spin frustration leads 
to a transition from the Neel state to a paramagnetic 
state at a critical on-site repulsion^r— . A metallic anti- 
ferromagnetic phase exists in a narrow region near this 
transition^—. Another possible ground state found in nu- 
merical studies^— due to spin frustration is the collinear 
antiferromagnetic (CAF) state in which the spin configu- 
ration is antiferromagnetic along one axis and ferromag- 
netic along the other axis. The transition from the Neel 
state to the CAF state in a Heisenberg model on a square 
lattice with both NN spin exchange J > and NNN spin 
exchange J' > was proposed at J' = and the 

exact phase diagram of this model is still under investi- 
gation. In the 2e?-Hubbard model, properties of the CAF 
state are largely unknown. 



In this work, we study properties of the CAF state in a 
half-filled Hubbard model on a square lattice at zero tem- 
perature. We find that there is a metal-insulator transi- 
tion in the CAF state. The magnetic moment displays a 
kink at this transition point. When U is small, the CAF 
state is metallic, and a van Hove singularity may be close 
to the Fermi energy When U is large, the CAF state is 
an insulator. The paramagnetic, antiferromagnetic, and 
CAF states are ground states in different regions of the 
phase space. There is a second-order transition between 
CAF and paramagnetic phases, and a first-order transi- 
tion between CAF and antiferromagnetic phases. Based 
on these results, a zero-temperature phase diagram is ob- 
tained, in comparison with previous studies^— . 



II. MEAN-FIELD THEORY OF THE CAF 
STATE 

The Hamiltonian of the Hubbard model with both NN 
and NNN hoppings is given by 

where Ci a and c| CT are electron annihilation and creation 
operators at site i with spin a, and — c^ ia Ci a is the 
number operator. The first two terms on r.-h.-s. of Eq. 
([1]) describe NN and NNN hoppings, and the last term 
describes the on-site repulsion. The square lattice has 
total N sites and the lattice constant is given by a. 

In the large U limit, the NN hopping leads to antifer- 
romagnetic spin exchange J = At 2 /U between NNs, and 
similarly J' = At' 2 /U between NNNs. The energy per site 
of the classical antiferromagnetic state is —(J — J')/2, 
while it is — J'/2 in the classical CAF state. Thus in- 
tuitively when J' > J/2 or cquivalently t' > t/y/2, the 
CAF state is preferred over the antiferromagnetic state. 
However, much is unknown when the repulsion U is of 
the same order as t or t', which will be investigated in 
this work. 

In both antiferromagnetic and CAF states, the square 
lattice can be divided into two sublatticcs. The order 
parameter is the same on each sublattice, but opposite 
on different sublattices. It can be generally written as 

(n it ) - (n^) = mcos(Q • n), (2) 

where is the coordinate of site i and to > is the 
magnetic moment. In the antiferromagnetic state, the 
wavevector Q is given by Q = (ir/a,ir/a); in the CAF 
state, Q = (0,-ir/a) or Q = (ir/a, 0). When m = 0, the 
order parameter is zero, and the system is in a paramag- 
netic state. 

In the mean-field approach, the on-site repulsion term 
in the Hamiltonian can be approximated by 

Un^Uix ~ U(n lt (riii) + n^n^) - {ni\){rii\)), (3) 

and the mean-filed Hamiltonian of the antiferromagnetic 
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or CAF states can be written as 

UN 2 , v^rV^/ Un x 
Umf= - — (n - m ) + 2^[)^( e k + —) c L c k* 

k a 

- (c kt c k+Qt - cJ, c k+Q I )] , (4) 



^-l c kt c k+Qt - c H c k+CUJJ> 

where e k = — 2t(cosk x a + cos k y a) + At' cos k x a cos k y a, 
n = (riif) + (n,4.) = 1 at half filling. 

The mean-field Hamiltonian can be diagonalized by a 
standard canonical transformation 



H 



MF 



UN 



n 2 — m 2 ) 



(5) 



where the k-summation is over the first Brillouin zone 
of a sublatticc, the quasi-particlc operators are given 
by a k t = WkCkf + v k c k+Q ^, 0-14 = u k c ki , - v k c k+ Qj., 
/?kt = -«kC k f + «kCk+Qt> /3kj. = «kCki + WkCk+Qi, and 
the coefficients are given by 
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FIG. 1: (a) Band structures of the CAF state at the metal- 
insulator transition point Umi = 4.996t for t' = 0.8t, (b) 
magnetic moment m versus repulsion U for t' = 0.8t. One 
kink of m appears at the transition point U = Umi, and the 
other kink appears at U = U v h = 3.74t where there is a van 
Hove singularity at the Fermi surface. 
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U|, = 1 — Vlr = — 



1 - 



£k - £k+Q 



V( £ k-ek +Q ) 2 +4A2 



with A = mU/2. The quasi-particles form two bands, 
with energies given by 



III. METAL-INSULATOR TRANSITION IN THE 
CAF STATE 

In a CAF state at half filling, c. g. with Q = (7r/a, 0), 
from Eq. ([6]) quasi-particle energies are given by 
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4 = 2 ( £k + ek +Q + C/n ) ± 2 V (6k " £k+q)2 + 4A ^ (6) 

The magnetic moment m can be determined self- 
consistently, 

m = 1 - £ k ) - % - 4)1, (7) 



where (i is the chemical potential. This self-consistency 
equation can be further written as, 

1 I flfc- g -)-g( M - £ +) 

- iv^ (8) 



k ^i(e k -e k+Q )2 + A 2 

Equation (jHJ can be solved together with the density 
equation 



= ]yrE'[^-4) + *(M-0] 



The total energy per site -E is given by 



(9) 



S = TcE' - 4) + MO" " ^)] - j(n 2 - ™ 2 )- 



(10) 



The self-consistency equation ((SJ) is equivalent to the 
energy-extreme condition dE/dm = 0. 



4 = iy^icos k x a — At' cos k x a cos k y a) 2 + A 2 

U 
+ 2 



2t cos k y a. 



(11) 



The minimum of the upper band is given by e^t- = — 2t + 
A + U/2 at k = (it /2a, 0), and the maximum of the lower 
band is given by e~ ax = 2i — A + J7/2 at k = (vr/2a, 7r/a). 
Therefore when i > A/2, the CAF state is metallic; when 
< < A/2, the system is an insulator with a band gap 2A — 
At. For t' = 0.8t, this metal-insulator transition takes 
place at a critical repulsion Umi = 4.996t, as shown from 
band structures in Fig. [Ha). When U > Umi, the CAF 
state is an insulator; when U < Umi, it is a metal. At the 
metal-insulator transition point, the magnetic moment m 
displays a kink as shown in Fig. [Ub). 

The other kink of m in Fig. QJb) is due to the van 
Hove singularity. From Eq. (jlip . when t'/t > 0.5, there 
is always a van Hove singularity in each band, located 
at (0, 4h) and (0, fc~ H ) in upper and lower bands respec- 
tively, kf H a = arccos[t/2i' ± t 2 & / (AtW At' 2 - t 2 )]. The 
density of states (DOS) diverges logarithmically at these 
points. In the insulating CAF state, the chemical poten- 
tial is always between the two bands at half filling, and 
van Hove singularities are not important. However, in 
the metallic CAF state, the chemical potential can reach 
the van Hove singularity in the upper band at certain 
repulsion U v yi, whereas the energy of the van Hove sin- 
gularity in the lower band is always less than chemical 
potential. 
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FIG. 2: Magnetic moment m versus repulsion [7 near C/ v h 
in the CAF state for (a) t'/t = 0.8 and (b) t'/t = 1. Total 
energies per site versus m axe in the insets, where Ep is the 
energy per site of the paramagnetic state. Band structures 
and DOS at U = U vH for (c) t'/t = 0.8 and (d) t'/t = 1. 



At E/vH, the magnetic moment m can exhibit either 
a kink as shown in Fig. [2] (a), or a jump as shown in 
Fig. [5] (b). For t'/t = 0.8, the magnetic moment m is 
continuous near £7 v h, because the total energy £ as a 
function of m has only one local minimum, as shown in 
the inset of Fig. [2(a). In contrast, for t'/t = 1, two local 
minima of E appear, as shown in the inset of Fig. [5] (b), 
and at C/ v h the magnetic moment is discontinuous. The 
van Hove singularity at t/ V H can be clearly seen from the 
band structure and the DOS divergence in Fig. H](c) and 
(d). In our numerical calculation, the discontinuity of m 
at U vli appears when t'/t > 0.84. When t'/t < 0.84, the 
magnetic moment m exhibits a kink at C/ v h- 



IV. PHASE DIAGRAM AT HALF FILLING 

The CAF state is not always energetically favored 
over antifcrromagnetic or paramagnetic states. It is well 
known that when there is only NN hopping, the ground 
state is the antifcrromagnetic state. When the NNN hop- 
ping is finite, electron spins are frustrated. In the strong 
coupling limit, U 3> t and U ^ i', in both CAF and anti- 
ferromagnetic states, the lower band is fully occupied and 
the upper band is empty From Eq. ((5]), the energy per 
site of CAF state is approximately E w -(2t 2 - it' 2 )/U, 
and in the antiferromagnetic state E « —At 2 /U. There- 
fore, when t'/t > l/y/2, the CAF state has lower en- 
ergy than the antifcrromagnetic state, consistent with 
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FIG. 3: (a) Critical interaction strength at the paramagnetic- 
antiferromagnetic transition as a function of t'/t. (b) The 
energy of the antiferromagnetic state 75n as a function of 
the magnetic moment ttin for different repulsion strength for 
t'/t = 0.05, with ttin = at the transition, (c) En for 
t'/t = 0.1, with m,N 7^ at the transition. The reference en- 
ergy Ep is the energy of the paramagnetic state for mN = 0. 



the J — J'— Heisenberg-modcl result in Ref.— . 

In the opposite limit, [f < f and U <§; t', the ground 
state is a metallic paramagnetic state. A continuous 
phase transition between paramagnetic and CAF states 
can take place at a critical repulsion Up which can be 
determined from Eq. (JSj) by setting A = 0, 



1 2 x - 

?7p ~ 



f e(jjL-e^)-9(ji-e+) 

kk - £ k+Q| 



(12) 



where = {U + e k + £k+Q ± |ek - £k+Q|)/2. The 
transition point between paramagnetic and antifcrromag- 
netic states can also be determined similarly. When 
t'/t < 0.707, the paramagnetic- antiferromagnetic phase 
transition takes place first with the increase of the re- 
pulsion U; when t'/t > 0.707, the paramagnetic- CAF 
transition occurs first. 

The paramagnctic-antifcrromagnetic transition was 
first found by Lin and HirscbA In later studies, it was rec- 
ognized that this transition can be either a first-order or 
second-order transition depending on the hopping ratio 
t' /1&z. A narrow region of a metallic antiferromagnetic 
state may exist near the transition^—. Our numerical 
results show that the second-order transition takes place 
not only in the region with t'/t > 0.38^, but also in the 
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FIG. 4: The mean-field phase diagram of the 2d-Hubbard 
model at half filling. AFI denotes antiferromagnetic in- 
sulator, AFM denotes antiferromagnetic metal, CAFI de- 
notes collinear antiferromagnetic insulator, CAFM denotes 
collinear antiferromagnetic metal, and PM denotes paramag- 
netic metal. 



region with t' ft < 0.08, as can be seen from the vanish- 
ing magnetic moment mN of the antiferromagnetic state 
at the transition in Fig. H When 0.08 < t'/t < 0.38, 
the magnetic moment mN is finite at the transition, indi- 
cating that it is a first-order phase transition. Near the 
second-order paramagnetic-antiferromagnetic transition, 
there is a tiny region of metallic antiferromagnetic state; 
near the first order transition, the antiferromagnetic state 
is always an insulator. 

By comparing energies of paramagnetic, antiferromag- 
netic, and CAF states, we have obtained the mean-field 
phase diagram of the 2d-Hubbard model at half filling as 
shown in Fig. [4] These three phases meet at a tricritical 
point (f = 0.707t,U = 4.01£). The paramagnetic state 
is always the ground state when the repulsion is small 
enough. When < t'/t < 0.707, there is a transition 
between paramagnetic and antiferromagnetic phases at 
a critical repulsion, which can be seen from the energy 
comparison shown in Fig. [5ja). When < t'/t < 0.08 or 
0.38 < t'/t < 0.707, it is a second order transition; when 
0.08 < t' /t < 0.38, it is a first-order transition, as shown 
in Fig. H 

In the region with 0.707 < t'/t < 0.774 and U > U P , 
CAF and antiferromagnetic states have close energies, 
and there is a reentrance from the CAF state to the 
antiferromagnetic state as can be seen from the energy 
comparison shown in Fig. Efb). There is a transition 
between the metallic CAF and antiferromagnetic states 
at a critical repulsion Un larger than Up. However the 
metallic antiferromagnetic region is very small. When the 
repulsion further increases, the antiferromagnetic state 
becomes insulating, and eventually a transition between 
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FIG. 5: Comparison between CAF-state energy Eq and 
antiferromagnetic-state energy _Z5n for: (a) t'/t — 0.65, (b) 
t'/t — 0.76 and (c) t'/t = 0.8, where Ep is the paramagnetic- 
state energy. 

insulating antiferromagnetic and CAF states takes place 
at another critical repulsion Uq- The transitions be- 
tween antiferromagnetic and CAF states at Uq and 
are the first-order transitions, because antiferromagnetic 
and CAF order parameters are discontinuous. Since the 
parameters U, t and t' are of the same order at these tran- 
sitions, fluctuations may modify these phase boundaries. 
In the phase diagram in RefJ^, the antiferromagnetic and 
CAF phases are separated by a very small paramagnetic 
region in between. In Refi^, the CAF-antifcrromagnctic 
transition was found near t'/t ss 0.78 without reentrance. 
The detail of this transition needs to be resolved in future 
studies. 

For t'/t > 0.774 and U > U P , the CAF state has al- 
ways lower energy than the antiferromagnetic state, as 
shown in Fig. [5tc). In the metallic CAF state, the mag- 
netic moment m has a kink for t'/t < 0.84 and a disconti- 
nuity for t'/t > 0.84, when the van Hove singularity is at 
the Fermi surface. The metal-insulator transition takes 
place inside the CAF state at repulsion Umi where also 
the magnetic moment m has a kink . 

V. DISCUSSION AND CONCLUSION 

The location of the CAF phase in the mean-field phase 
diagram is in qualitative agreement with those obtained 



5 



in other approaches^— . We would like to stress that the 
mean-field method is more accurate when the on-site re- 
pulsion is relatively small. It is inadequate for treating 
spin-liquid or superconducting states ^ 10 ' 16 due to strong 
correlation. Spin frustration can significantly change 
the magnetic excitations in the antiferromagnetic stated. 
Beyond half filling, ferromagnetic state, incommensurate 
spin-density-wave state, and phase separation^ between 
these states can also appear. 

In conclusion, we studied the CAF state in a half-filled 
2d Hubbard model with NN and NNN hoppings at zero 
temperature. We found that a metal-insulator transition 
takes place inside the CAF phase at a critical on-site re- 
pulsion Umi- In the metallic CAF state, there is a kink or 
discontinuity in the magnetic moment when the van Hove 



singularity is at the Fermi surface. At zero temperature, 
the CAF, antiferromagnetic, and paramagnetic states 
meet at a tricritical point (t' = 0.707t, U = 4.01t). There 
is a first-order CAF-antiferromagnetic phase transition, 
and a second-order CAF-paramagnetic phase transition. 
The mean-field phase diagram of the 2d Hubbard model 
is obtained. 
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